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Abstract 

We discuss the relationship between the RR interval spectral analysis and 
the spectral analysis of the corresponding ECG signal from which the RR 
intervals were evaluated. The ECG signal spectrum is bounded below the 
frequency ffi by using an electronic filter and sampled at rate larger than 2fB, 
thus excluding aliasing from spectral analysis. A similar procedure cannot 
be applied to the RR interval spectral analysis, and in this case aliasing is 
possible. One of our main effort in this paper is devoted to the problem of 
how to detect aliasing in the heart rate spectral analysis. In order to get an 
insight we performed an experiment with an adult man, in which the ECG 
signal was detected in a case where the breathing rate was larger than half the 
heart rate. A constant breathing rate for time intervals exceeding 5 minutes 
was monitored with good accuracy using a special breathing procedure. The 
results show distinctively a very sharp peak in the spectral analysis of the 
ECG signal and corresponding (diffused) aliasing peaks in the RR interval 
spectral analysis. 

New method of dealing with unevenly sampled data was developed which 
has interesting anti-aliasing properties. There are indications that the VLF 
peaks of the RR spectrum are originated by aliasing. Some of the LF peaks 
may have the same property. 
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I. INTRODUCTION 

The R-R interval spectral analysis is usually based on heart rate data collected in two 
ways. In one method the data are collected by analog to digital conversion of the ECG 
signal and computer evaluation of the R-R intervals from the ECG signal. In the second 
method, devices are used whose output is the R-R interval alone. The advantage of the first 
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method is the control of accuracy and flexibihty of the evaluations. The second method has 
the advantage of storing smaller amount of data, and it can be easily used on-line. 

In the first method, usually the number of collected data (sampled ECG signal) is 
of two to three orders of magnitude larger than the R-R interval data. Thus if only R-R 
interval is analyzed a large amount of data is unused. In this paper we are trying to take 
advantage of the ECG sampled signal and to derive new information in addition to the 
conventional R-R interval analysis lH,, 0, [|], [Q. 

The ECG signal spectrum is bounded below the frequency fs by using an electronic 
filter and sampled at rate larger than 2fB, thus excluding aliasing from spectral analysis. 
A similar procedure cannot be applied to the R-R interval spectral analysis, and in this case 
an aliasing is possible. One of our main efforts in this paper is devoted to the problem of 
how to detect aliasing in the R-R interval spectral analysis. 

In order to get an insight, we performed an experiment, in which the ECG signal 
of one of the authors (A.G) was detected while the breathing rate was larger than half the 
heart rate. A constant breathing rate for a time exceeding 5 minutes was monitored with 
good accuracy using a special breathing procedure with a metronome. The results show 
distinctively a very sharp peak in the spectral analysis of the ECG signal and corresponding 
(diffused) aliasing peaks in the R-R interval spectral analysis. 

The spectral analysis of the ECG signal was performed with the standard FFT proce- 
dures. The spectral analysis of the R-R intervals was performed with several techniques in 
order to take into consideration that the data were unevenly sampled. This is presented in 
section 2. In section 3 we discuss the possibility of aliasing in the spectral analysis of the 
R-R intervals. In section 4 we compare power estimations of ECG's and R-R intervals of 3 
experiments. In section 5 we analyze the results. In section 6 summary and conclusions are 
presented. 
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II. SPECTRAL ANALYSIS OF UNEVENLY SAMPLED DATA 



The methods of spectral analysis are well developed for evenly sampled data [§[,[§[• The 
R-R interval data are unevenly sampled in time. In most cases an analysis is performed 
with respect to beat numbers which are evenly spaced. We will below justify this method 
using least squares principles. But as was recently indicated by Laguna et Al. [[7|, the 
resampling of data is causing the appearance of additional harmonics. They recommend to 
use a method developed by N. R. Lomb [§] . The errors of resampling the beats can, to large 
extent, be overcome by using a cubic spline interpolation. In this work we are suggesting 
a new method of treating unevenly sampled data, which, unexpectedly, gave good results 
beyond the Nyquist frequency. 

A. Analysis according to beat numbers 

Let us assume that the RR intervals are given at unevenly sampled times tn, with the 
values s , where n is the beat number, n = 1- ■ ■ N. Let us divide the interval [ti,tAr] 
into equal subintervals 

and let us generate in the interval [ti,tAr] evenly sampled times: 

r„ = (n- l)Ar + ti. (2) 

We will use the discrete time Fourier transform (DFT) for a basis formed from the evenly 
sampled times r„. We will assume that 

1 ^ 

s{tn) = TT^^Sk exp (iUkTn) , c^fc = 27r (A; - 1) / (iVAr) . (3) 
k=i 

The coefficients Sk will be determined by minimizing the expression 



N 



n=l 



1 ^ 

S(tn) - ir^'^Sk exp {iUkTn 

k=l 



1 ^ 

S(tn) - Jt'^SI exp {-iUkTn 
k=l 



(4) 



with the result 

N 

Sk = ^s (t„) exp {-iuJkTn) . (5) 

n=l 

Eqns. ^ and ^ can be handled easily with standard FFT programs. This is the usual 
procedure which is adopted in most of the papers dealing with R-R interval analysis. ^ ^ 

B. Other Methods 

FFT can be applied more efficiently if the unevenly sampled data are interpolated at 
evenly spaced intervals of Eq. ||. The cubic spline interpolation is one of the good ways to 
do it. 

The Lomb method was extensively analyzed in ref . [0 . We give here only the formulae 
in the form of the Lomb normalized periodogram 



Px M 



En=l (tn) - S] COS (cUfc (t„ - t)) E„=1 [s [tn) " s] siu (Wfc {t^ - t)) 



+ 



2 ^ 



1 En=l C0S2 {Uk {tn - t)) siu^ {Uf, {tn ' t)) 

(6) 



where s and a are the mean and variance of the data and the value of r is defined as 

tan(2c.,r) = ^r-^^^4^ (7) 

E^=iCos(2a;,t„) ^ ^ 



C. Non-Uniform Discrete Fourier Transform (NUDFT) 

We present here a new method of treating unevenly spaced events which we call the 
"non-uniform discrete Fourier transform" (NUDFT). 

Let us assume that s (r„) are the exact values of the signal at the points given by Eq. |^. 
The corresponding DFT is 

N 

Sk='^S {Tn) exp {-iuJkTn) ■ (8) 
n=l 
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Our aim is to find a good approximation to this expression in terms of tlie unevenly 
sampled signal 

We start with the Euler summation formula 

N 



1 f^N 1 At 

E/(^n) = ^/^ /(r)rfr+-[/(n) + /M] + — [/'M-/'(n)] + 0(Ar2) (g) 
and make the following decomposition of the integral on the right hand side of Eq.^: 

/ f{T)dr= f{T)dT+ f(r)dT+---+ f{T)dT (10) 

and approximate each of the integrals on the right hand side with the trapezoidal rule: 
/J / (r) dT = l[f (ti) + / ih)] (t2 - ti) + ■ ■ ■ + ^ [/ (tN-i) + f (tiv)] {tN - tN-i) + O (At) 

(11) 

From Eqs. |^ and |TT] we obtain: 

E / = ^ {[/ (^l) + / (^2)] (t2 - tl) + ■ ■ ■ + [/ + / (t^)] {t^ - t^-l)} 

+ l[f{h) + fM]+0{Ar). (12) 



When the t„ are equally spaced Eq. |T2] becomes an identity with the O (Ar) = 0, 



therefore it seems to us that Eq. |T2| is satisfied with an higher accuracy than just O (Ar) . 
Eq. can be applied to approximate Eq. ^ with the substitution 

/ (tn) = S (tn) exp {-iUktn) , (13) 

and the final result, the approximation to Eq. ^, after rearranging the terms, becomes: 

N 

Sk = Y^ CnS (tn) exp {-iUktn) + O (At) , (14) 

n=l 

where 
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with the inverse formula 



1 ^ 



_ Ar+f2-tl 
•"1 2Ar ' 

•"2 2Ar ' 



tN—tN-2 
CN-1 - 2Ar ' 

„ Ar +f Af — tjv-i 

- 2A7 ' 



(15) 



= ]^ E "^fc exp (zc^fcr„) + O (Ar) , c^^ = 27r (A; - 1) / (NAr) , (16) 
fc=i 

which is an interpolation formula for s (tn) at the evenly spaced points ri ■ ■ ■ ttv- 



III. ALIASING 



Aliasing is a result of undersampling and is a well known phenomenon. In ref. ^ aliasing 
was looked upon from the point of view of symmetry. It is an example of wrong symmetry, 
and as such should be given more attention. It is the outcome of an incomplete basis. It 
was found in ref. , that for evenly sampled data with a sampling rate fs, the spectral 
amplitude S{f) evaluated with FFT, has the following symmetry properties 

\S (/)| = \S{f± fs)\ = \S (-/ ± fs)\ = \S (±/ ± nfs)\ , (17) 

where / is the frequency and n is an arbitrary integer. 

In order to avoid the aliasing symmetry of Eq. [l^, the frequencies should be bounded 
by the Nyquist frequency (denoted here by fs) according to 

fB = y. (18) 

The ECG signal was sampled with sampling rate 250 Hz, and an electronic filter was 
applied, which have eliminated practically all frequencies above 32 Hz, thus aliasing can not 
occur at frequencies below 125 Hz or even below 32 Hz. The R-R intervals were calculated 
directly from the ECG signal. The sampling rate for R-R intervals can be defined only for 
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evenly sampled data, for the methods which interpolates the unevenly sampled data, or one 
can consider the average sampling rate from Eq. |I|, in both cases 

fs = 1/At = 2fM, (19) 

where Jn is the Nyquist frequency for the R-R intervals. As the EGG signal contains 
frequencies much grater than /jv, and the R-R intervals are derived from the EGG signal, 
one can not be sure that the spectral analysis of the R-R intervals is free from aliasing. As a 



matter of fact there are indications of aliasing in some rare cases. ||T0[ ||TT| |[T2| ||T3[ |Tj] [15 
One way to identify aliasing is to change the sampling rate and follow the changes in the 
spectrum. Unfortunately, for the R-R intervals, one can not speak about a definite sampling 
rate, but rather can consider a distribution of sampling rates. The changes in sampling rate 
required to observe aliasing are of the same order as the fluctuations in the sampling rate. 
Therefore in practice it is almost impossible to observe consistent changes in the spectrum 
slightly changing the heart rate. 

Other possibility of detecting aliasing is by comparing the heart rate spectrum with 
the EGG signal spectrum. Marked differences below the Nyquist frequency for the power 
distribution of the RR intervals compared to the EGG signal power distribution in the same 
range may indicate aliasing. But we do not have yet a sound basis to treat this problem. 

We have devised an experiment which definitely demonstrates the aliasing in the R- 
R intervals spectrum. To the best of our knowledge this is the first experiment in which 
one can exactly know the correct frequency above the Nyquist frequency and can follow 
the development of the aliasing, which appears to be diffused to great extent because the 
symmetry of Eq. |1^ is represented not by one sampling rate but by a distribution of sampling 
rates, as the R-R interval is unevenly sampled. 

Below we describe 3 experiments. One of them was devised to demonstrate aliasing and 
the other two for learning about the relations between the R-R interval spectrum and the 
spectrum of the EGG signal. 



8 



IV. THREE EXPERIMENTS 



We present below results of three experiments. In the first experiment the ECG signal 
was collected in a normal resting state. The aim of this experiment was to compare the ECG 
spectrum with the R-R intervals spectrum. In the second experiment very slow breathing 
was monitored at a rate of 0.04 Hz. Again the ECG and R-R interval spectra were compared. 
In the third experiment very fast breathing was accurately monitored at the rate of 74/min 
and 84/min. These respiratory rates were above half of the heart rates thus allowing to 
observe in detail the development of aliasing. 

A. The First Experiment 

In this experiment (linked with the names of Zahi and Ori, where the second is one of 
the authors: O.G) which was done in normal, resting conditions, we compare the power es- 
timation of the R-R interval and the ECG signal, from which the R-R interval was obtained. 
The ECG signal was sampled at a rate of 250 Hz. Stable intervals of 7 minutes duration 
were chosen for analysis. 

In Fig. la the power distribution of the ECG signal of Zahi is depicted. The attenuation 
of the power with increasing frequency above 12 Hz is due to the action of an electronic 
filter. Above 32 Hz the contribution is practically zero. The average heart rate was 0.97 Hz. 
The above results were zoomed to the interval [0-12] Hz in Fig. lb. One can see distinctively 
the peak around the average heart rate and the higher harmonics of this peak. The second 
harmonic is missing, but the third, fourth, fifth and sixth are distinctively visible, higher 
harmonics became more and more smeared and indistinguishable above the sixth harmonic. 
One should also note the large difference in power in the heart rate range, below the Nyquist 
frequency of 0.49Hz, which is much smaller compared to the peak around the average heart 
rate (0.97 Hz). 

The power distribution of the RR intervals in the range [0-0.5] Hz was computed ac- 
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cording to the methods discussed in section 2 and are presented in Figs. 2a (DFT, beat 
number analysis), 2b (Spline interpolation), 2c (NUDFT). For comparison also the power 
distribution of the ECG signal in the above range is presented in Fig. 2d. 

The results of Figs. 2a, 2b, and 2c are quite similar, but the sphne interpolation (Fig. 
2b) and the NUDFT (Fig. 2c) are practically identical. The three graphs show the structure 
commonly found in the power estimation analysis of RR intervals, namely the existence of 
the "high frequency" (HF), "low frequency" (LF) and the "very low frequency" (VLF) 
peaks. The ECG spectrum shows qualitatively the same structure (but not a quantitative 
agreement), except that the ECG spectrum is highly suppressed below 0.04 Hz, in the VLF 
region, indicating a possibility of aliasing in this region in the RR analysis. 

In Figs. 3 and 4a-4d the results of Ori are presented. The conclusions are similar to those 
of Zahi, except that in the ECG spectrum both VLF and LF peaks are missing, indicating 
the possibility of ahasing in these regions for the RR analysis. Also in the ECG spectrum 
of Ofek, Fig. 5b the VLF and LF, present in Fig. 5a, are missing. VLF is missing in J.C.'s 
ECG spectrum (see Figs. 6a-6b). 

B. The Second Experiment 

In this experiment (linked again with the name Ori) we have checked the ECG spectrum 
near the VLF region, as the VLF was absent in the ECG spectrum for the resting state in 
the first experiment. The question was whether such a result persists in all ECG spectra. 
Therefore we have probed the VLF region by monitoring very prolonged breathing with a 
rate of 0.04 Hz. For the spectrum of RR intervals we found that the DFT, Spline interpo- 
lation and NUDFT gives similar results, and again NUDFT was practically identical to the 
spline interpolation. Therefore we present only the results of NUDFT, which are presented 
in Fig. 7a. For comparison the spectrum of the ECG signal is given in Fig. 7b. In Fig. 7a 
one can see a very clean pattern of a peak at 0.04 Hz and its higher harmonics. In Fig. 7b 
one can see a similar but somewhat diffused pattern. Thus this experiment indicates that 
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similar respiratory patterns exists in both the RR as well as in the EGG signal. 

C. The Third Experiment 

In this experiment (linked to the name Alex, who is one of the authors: A.G) very fast 
breathing was accurately monitored at the rate of 74/min and 84/min respectively. These 
rates were well above half of the average heart rate thus allowing to observe in detail the 
development of aliasing. In Fig. 8 the EGG spectrum is dominated by the very high and 
narrow peak at the frequency /i = 1.23AHz, also its higher harmonics can be distinctively 
seen. The frequency /i is just the breathing frequency 74/min. In the same figure one can 
also see the diffused peaks near the average heart rate frequency of 1.636 Hz and its higher 
harmonics. One should observe aliasing at about 1.636— /i = 0A02Hz. Indeed one can 
see diffused peaks around that frequency in Fig. 9a, which displays the power estimation of 
the RR intervals using the NUDFT (which below the Nyquist rate is similar to the spline 
interpolation). The width of this region can be estimated by noting that the RR intervals 
have different instantaneous sampling rates which are equal to the inverse of the RR interval 
time. In Fig. 10 we have calculated the distribution of the sampling rates by dividing the 
frequency region into 100 beans. We have shifted that distribution by subtracting fi. As 
one can see the results are confined approximately to the region 0.32-0.47 Hz. Indeed the 
aliasing peaks of Fig. 9a appear in this region. The pictures below the Nyquist frequency 
are very similar for the DFT, NUDFT, the spline interpolation and the Lomb method (Fig. 
9b) with a similar aliasing behavior. 

In principle the NUDFT and the Lomb methods should not be used above the Nyquist 
frequency. Surprisingly enough we have found that both methods have a sharp peak at /i, 
as can be seen in Figs. 9a and 9b. Both methods do not have the aliasing symmetry of the 
DFT as given by Eq. therefore the results are not symmetric with respect to the Nyquist 
frequency (half the sampling rate), as it is satisfied, for example, in the case of the spline 
interpolation. We have found an exact result at /i and a diffused ahasing around 0.4 Hz. It 
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is interesting to note that both methods give almost the same result below and above the 
Nyquist frequency. One can interpret the appearance of the sharp peak at /i as a result of 
a partial destruction of aliasing symmetry due to uneven samplings. 

Similar results for the breathing frequency 84/min are presented in Figs. 11-12. 

V. SUMMARY AND CONCLUSIONS 

The EGG signal spectrum is bounded below the Nyquist frequency fs by using an 
electronic filter and sampled at rate larger than 2iB, thus excluding aliasing from spectral 
analysis. A similar procedure cannot be applied to the RR interval spectral analysis, and 
in this case an aliasing is possible. One of our main efforts in this paper was devoted to the 
problem of how to detect aliasing in the R-R interval spectral analysis. 

In order to get insight into this problem three experiments have been analyzed. In the 
first experiment the EGG signal was collected in a normal resting state. The aim of this 
experiment was to compare the EGG spectrum with the R-R interval spectrum. In the 
second experiment very slow breathing was monitored at a rate of 0.04 Hz. Again the EGG 
and R-R interval spectra were compared. In the third experiment very fast breathing was 
accurately monitored at the rate of 74/min and 84/min respectively. These respiratory 
rates were above half of the heart rates thus allowing to observe in detail the development 
of aliasing. 

The experiments which were described above led us to the following conclusions: 

1. The spectral analysis of the EGG signal is more sensitive and accurate compared to 
the R-R interval spectral analysis and is free from ahasing. Still in the present stage 
it contains too much information to be of practical use. Efforts should be made to 
understand what will be the best way to extract information (not related to the heart 
condition alone as in the standard analysis of EGG) about the external influences on 
the heart signal. 
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2. We have conducted an experiment which gave a clear insight about the mechanism of 
ahasing in the R-R interval spectrum. The very sharp peak in the spectrum of the 
ECG signal, which came as the result of enforced quick breathing, reappeared as a 
diffused signal in the RR spectrum. The extension of the diffuseness agrees with the 
extension of the sampling rates of unevenly sampled data.. 

3. The VLF peak observed in the R-R interval spectrum is usually missing in the ECG 
spectrum. This lead us to suspect that the VLF observed in the RR spectrum has its 
origin in aliasing. 

4. In some cases the LP peak does not show up in the ECG spectrum. This led us to 
suspect that part of the LF peak is of aliasing origin. 

5. Unlike in electronic devices, it is very difficult to devise procedures to detect aliasing in 
humans. In electronic devices aliasing can be easily detected by changing the sampling 
rate. In humans the fluctuations of the heart rate are of the same order as the required 
changes in the sampling rates. It will be an important task to develop a proper 
procedure for detecting aliasing in humans. 

6. We have developed a new technique for spectral analysis for unevenly sampled data 
called non-uniform discrete Fourier transform (NUDFT). When employed to the RR 
data, below the Nyquist frequency, it gave similar results as those obtained by interpo- 
lating the data with a cubic spline. Above the Nyquist frequency, the correct peak in 
the spectrum was detected with great accuracy. A similar result was obtained with the 
recently rediscovered Lomb method. We interpret this unexpected result by a partial 
destruction of aliasing symmetry in both methods. More efforts should be made in 
order to understand the anti-aliasing properties of the above methods. 

7. We consider aliasing to be a wrong symmetry, resulting from the use of an incomplete 
basis, which has intrinsic symmetries inconsistent with the properties of the signal. 
Aliasing can be partially removed by reducing the symmetry of the basis. 
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Figure Captions 

Figure 1. The relative power of the ECG signal of Zahi, a) in the spectral range of 
0-36 Hz, b) in the spectral range of 0-12 Hz. 

Figure 2. The relative power computed (from the ECG signal of Zahi) by four different 
methods, in the spectral range of 0-0.5 Hz, a) by DFT, b) by spline interpolation of 
the RR data, by NUDFT, d) from the ECG signal. 

Figure 3. The relative power of the ECG signal of Ori. 

Figure 4. The relative power computed (from the ECG signal of Ori) by four different 
methods, in the spectral range of 0-0.52 Hz. a) by DFT, b) by sphne interpolation of 
the RR data, c) by NUDFT, d) from the ECG signal. 

Figure 5. The relative power computed (from the ECG signal of Ofek) by two different 
methods, in the spectral range of 0-0.6 Hz, a) by spline interpolation of the RR data, 
b) from the ECG signal. 

Figure 6. The relative power computed (from the ECG signal of J.C.) by two different 
methods, in the spectral range of 0-0.46 Hz, a) by sphne interpolation of the RR data, 
b) from the ECG signal. 

Figure 7. The relative power computed (from the ECG signal of Ori with breathing 
rate of 0.04 Hz) by two different methods, in the spectral range of 0-0.62 Hz, a) by 
NUDFT, b) from the ECG signal. 

Figure 8. The relative power of the ECG signal of Alex with a breathing rate of 1.234 
Hz. 

Figure 9. The relative power computed (from the ECG signal of Alex with a breathing 

rate of 1.234 Hz) by two different methods, in the spectral range of 0-1.5 Hz, a) by 
NUDFT, b) from the ECG signal. 
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• Figure 10. A 100 bin histogram of the heart rates of Alex which are subtracted by the 
breathing rate of 1.234 Hz. 

• Figure 11. The relative power of the ECG signal of Alex with a breathing rate of 1.404 
Hz. 

• Figure 12. The relative power computed (from the ECG signal of Alex with a breathing 
rate of 1.404 Hz) by two different methods, in the spectral range of 0-1.6 Hz, a) by 
NUDFT, b) from the ECG signal. 
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